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Abstract Tunneling spectroscopy at surfaces of unconventional superconductors 
has proven an invaluable tool for obtaining information about the pairing symmetry. 
It is known that mid gap Andreev bound states manifest itself as a zero bias conduc- 
tance peak in tunneling spectroscopy. The zero bias conductance peak is a signature 
for a non-trivial pair potential that exhibits different signs on different regions of 
the Fermi surface. Here, we review recent theoretical results on the spectrum of 
Andreev bound states near interfaces and surfaces in non-centrosymmetric super- 
conductors. We introduce a theoretical scheme to calculate the energy spectrum of 
a non-centrosymmetric superconductor. Then, we discuss the interplay between the 
spin orbit vector field on the Fermi surface and the order parameter symmetry. The 
Andreev states carry a spin supercurrent and represent a helical edge mode along the 
interface. We study the topological nature of the resulting edge currents. If the triplet 
component of the order parameter dominates, then the helical edge mode exists. If, 
on the other hand, the singlet component dominates, the helical edge mode is absent. 
A quantum phase transition occurs for equal spin singlet and triplet order parameter 
components. We discuss the tunneling conductance and the Andreev point contact 
conductance between a normal metal and a non-centrosymmetric superconductor. 
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In this chapter, we will discuss the surface and interface properties of non-centrosym- 
metric superconductors 1 35 1 focusing on the tunneling conductance. Since the early 
sixties tunneling spectroscopy has played an important role in gathering informa- 
tion about the gap function of conventional superconductors 0J. I n the context of 
unconventional superconductivity tunneling spectroscopy appeared as an important 
tool to probe the internal phase structure of the Cooper pair wave functions EJGl • 
Surface states with sub-gap energy, known as Andreev bound states (ABS) [|4]-[7] 
provide channels for resonant tunneling leading to so-called zero-bias anomalies in 
d//dV. Zero-bias anomalies observed in high-temperature superconductors showed 
the presence of zero-energy bound states at the surface, giving strong evidence for 
c/-wave pairing [2-6]. Similarly the tunneling spectrum observed in Sr2RuC»4 is con- 
sistent with the existence of chiral surface states as expected for a chiral p-wave su- 
perconductor ll8l- [T0l . Zero bias conductance peaks due to Andreev bound states have 
been observed in numerous experiments, e.g. in high-r c cuprates 0TJ, Sr2Ru04 
I112II13I . UBen [14], CeCoIns [15|, the two dimensional organic superconductor K- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br [16] and PrOs 4 Sbi 2 ED- Andreev bound states have 
also been observed in the Balian-Werthammer phase of superfluid 3 He [18]. The 
study of Andreev bound states in unconventional superconductors and superfluids 
has emerged as an important phase sensitive probe. 

In section |2] we present the theory for Andreev spectroscopy using Bogoliubov 
wave function technique in Andreev approximation. Starting with superconductors 
exhibiting li-wave or /j-wave pairing, we proceed with non-centrosymmetric su- 
perconductors. In section|3] we develop the theoretical tools for describing Andreev 
spectroscopy in non-centrosymmetric superconductors in the framework of Nambu- 
Gor'kov Green's functions within the quasiclassical theory of superconductivity. 



2 Andreev spectroscopy in unconventional superconductors 

2.1 Andreev conductance in s- and d-wave superconductors 

We discuss first the example of zero-bias resonant states at the interface of a normal 
metal/spin-singlet d- wave superconductor junction. In general, the pair potential can 
be expressed in terms of two coordinates, x and x', as A(x.x' ). In uniform systems it 
only depends on the relative coordinate x — x', and a Fourier transform with respect 
to it yields A (k) with relative momentum k. For illustrative purposes, we assume in 
the following a cylindrical Fermi surface and concentrate on two-dimensional sys- 
tems. The pair potential for spin-singlet cZ-wave pairing is A (9) = Aqcos(29), with 
e l9 = (kx + iky)/ \k\, while the corresponding spin-singlet s-wave one is isotropic, 
A(6) = Aq. The bulk quasiparticle density of states normalized by its value in the 
normal state is given by 
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Pb(E) = - J* &6 Po {E,A{0)), po(E,A(6)) = E =. (1) 

Kjo ^£ 2 -4 o 2 cos 2 (20) 

For a spin-singlet c/-wave superconductor this quantity behaves linearly at low ener- 
gies, p (E ) oc \E | . As shown below if the angle between the interface normal and the 
lobe direction of the li-wave pair potential has a nonzero value a with < a < n/2, 
then the resulting tunneling conductance <Jj(E) has a zero bias conductance peak. 

The Andreev conductance for a normal metal/insulator/spin singlet s-wave super- 
conductor junction is described by the model of Blonder, Tinkham, and Klapwijk 
(BTK) 1 19|. Within this model, Oj(E) at zero temperature is given by 

ME) « I fl (l + | a(E, 9) | 2 - | b(E, 9) | 2 ) (2) 

where a(E,9), and b(E,9) are probability amplitude coefficients for Andreev re- 
flection and for normal reflection, respectively. We apply the BTK model in the 
following to the case of spin-singlet d-wave pairing. 

We assume that the Fermi energy Ef is much larger than | A (9) , such that the Andreev ap- 
proximation can be applied to the Bogoliubov wave functions. For simplicity, we also assume 
equal effective masses and Fermi momenta in the normal metal and in the superconductor. The 
spatial dependence of the pair potential is chosen to be A(6)0(x) (with the Heaviside step function 
®). The insulating barrier at the atomically clean interface is modeled by a 5-function potential, 
V(x) = HS(x). Since the momentum parallel to the interface is conserved, the two component 
Bogoliubov wave function is given in Andreev approximation by 

n9 ' x) = ( V+(m) ex P(^cose)+ ( "-[^j) exp(-%xcos6) (3) 
where Uj(6,x) and Vj(6,x) with (j = +, — ) obey the Andreev equations 

Euj(0,x) = - [ - ajk * cos 9 j. _ HS(x)] u j (8,x)+A(e j )0{x)vj(8,x), 
L m ax J 

Evj(9,x) = [ - g/ ^ COs9 A_ g g(^)j v j (B,x) + A*(e j )0(x)u j {B,x), (4) 

with <j + = 1, 9 + = 0, and cr_ = — 1, 6 =n — 6. For a d-wave superconductor the corresponding 
effective pair potentials A(6±) are given by 

A(8+) = 4 cos(29-2a), 4(6_) = 4 cos(29 + 2a), (5) 

where the angle between the interface normal and the lobe direction of the rf-wave pair is a. The 
wave functions u±(6,x) and v±(6,x) resulting from Eqs. @ are obtained from by ansatz 



u + (6,x) 



q j exp(i5x) + a(E, 0) ( *j* J exp(— iSx) x < 



- {6 > X) J \c(EB)( V(E + n + )/2E \ (i , X>Q 

6) \ exp(-i^) y/(E - n + )/2E J eXp(l7+X) X> ° 

„_,e., M iK*,*)(;)e*p(-iS*) x<o 

v_(0,ac) J | ./„ m ?exp(i0_W(£-.Q-)/2£\ ,. . _ 



(6) 



(7) 
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Fig. 1 (Color online) Left: Tunneling conductance for an .v-wave superconductor, a: Z = 0, b: Z = 1 
and c: Z = 5. Right: Tunneling conductance for a (/-wave superconductor for Z = 5. a: a = 0, b: 
a = 0. 1257T and c: a = 0.257T. 



where we used the abbreviations 



^ ft = ^~ a , n ± = JE1-Al(6 ±) , exp W± ) = -^L. (8) 



hkf cos 9 hkf cos 6 V 4 ( 

With the help of appropriate boundary conditions, 

W(0,x) |, =0 _=f(9,*) \ x=0+ 

^(e,x)\ -^-w(e,x)\ =^ , p(9,x)| (9) 

dx v '\x=o+ dx y y Lv=o_ h 2 \x=o+ 
we obtain a(E, 9), b(E, 9), c(E, 9), and d(E, 6). 

The resulting conductance is 

o T (E)= (^J* * d6 D(6) Gr(E,6) cosflj /^^d8D(9) cos0 



i+D(e)|r + | 2 -fl(e)|r + n 
i-/?(0)r+r_ex P [i(0_-0 + )] | 2 

withf± = (E — £2±)/ | 4(0±) |. The quantities and/?(0) above are given by 



D(0) =4cos 2 0/(4cos 2 0+Z 2 ), /?(0) = 1 -£>(0), 

with injection angle and Z = 2mH /t^k? J2J. Choosing A(0±) = 4o reproduces 
the BTK formula for an s-wave superconductor. Typical line shapes of Gj(eV) with 
<?V = £ for s-wave and li-wave superconductors are shown in Fig. [T] The li-wave 
case is shown in Fig. Q2b)- As can be seen there, if the angle a deviates from 0, 
the resulting d//dV has a zero bias conduction peak (ZBCP) (curves b and c); the 
only exceptional case is a = 0, as shown in curve a. The width of the ZBCP is 
proportional to D, while its height is proportional to the inverse of D. The origin of 
this peak are mid gap Andreev bound states (MABS). The condition of the formation 
for Andreev bound states at the surface of an isolated cZ-wave superconductor (D — ► 
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Fig. 2 (a)+(c): order param- 
eter amplitude and (b)+(d): 
local density of states at the 
interface for a layered d- 
wave-superconductor/normal- 
metal junction. (a)+(b): a = 0, 
(c)+(d): a = k/4. The inter- 
face is at x = 0. The curves 
are for the indicated trans- 
mission coefficients Do. The 
temperature is T = 0.3r c , and 
the mean free path £ = lOijo- 
After Ref. tM. 
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0) is expressed by 

i = r+r_ex P [i(4>_-4> + )]. (ii) 

At zero energy = — 1 is satisfied, and consequently a MABS appears provided 
exp[i(0 + — 0_)] = — 1. For this case, on the superconducting side of the interface, 
the injected electron and the reflected hole experience a different sign of the pair 
potential. For a = n/4, there is a MABS independent of the injection angle. In this 
case, the energy dispersion of the resulting ABS, E^, is given by 

E h =0. (12) 

Finally, we comment on the effects of order parameter suppression near an sur- 
face or interface in a d-wave superconductor. In Fig.|2]we reproduce a self-consistent 
solution for a layered d-wave superconductor, showing that a strong order parameter 
suppression is always present for a = 0.25, whereas for a = in the tunneling limit 
the order parameter suppression can be neglected. The corresponding local density 
of states at the surface is shown in Fig.|2b) and (d). The interface is modeled by a 
5-potential as above, with a transmission D(9) = A) cos 2 0/(1 — Do sin 2 0), and the 
parameter Do is related to Z via Do = 1/[1 + (Z/2) 2 ]. 



2.2 Andreev conductance in chiral p-wave superconductor 

In this section, we discuss the tunneling conductance of a normal metal/chiral p- 
wave superconductor junction. There is evidence supporting the realization of spin- 
triplet pairing with broken time reversal symmetry in the superconducting state 
of S1-2RUO4 ll2TH26l . A possible symmetry is given by two-dimensional chiral p- 
wave pairing, where the pair potentials are given by At- * = A\\ = 0, At- j, = + = 
4o ex p(if}). In the following, 6 is measured from the interface normal. In the actual 
sample, the presence of chirality may produce chiral domain structures. A recent 
experiment is consistent with the presence of chiral domains [27 1 . Also, there are 
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several theoretical proposals to detect chiral domain structures [28 1. Here, for sim- 
plicity, we consider a single domain chiral p-wave superconductor. 

Since the z-component of the Cooper pair spin is zero, we can also use Eq. ( fTOb 
to obtain the tunneling conductance for normal metal/chiral p-wave superconductor 
junctions. Before discussing the tunneling conductance, we first consider the bulk 
local density of states (LDOS) of a chiral p-wave superconductor. In contrast to the 
spin-singlet d-wave pairing case, po(E,A(Q)) in Eq. (T]) is given by 

Po(E,A(6))=E/tJe2-AI 

It has a fully gapped density of state like in the spin-singlet s-wave case. 

We now discuss the condition when an ABS is formed at the surface of an isolated 
chiral p-wave superconductor. The bound state condition is given by J8][9] 

E+^/e 2 -A$ = -(e- yj& - A%\ rap(-2ie), (13) 

showing that the bound state level E\, satisfies 

E b {d)=A sm9 (14) 

Note that the ABS has a dispersion different from that in the c/-wave case with 
a = n/A. The presence of the edge state with a dispersion induces a spontaneous 
dissipationless current. 

As in the previous section we consider the tunneling conductance (Jj(E) in a 
normal metal/chiral p-wave superconductor junction, which is shown in Fig. |5J a )- 
As can be seen, for Z = 0, the line shape of conductance is identical to that of a 
spin-singlet s-wave superconductor (see curve a), whereas with increasing Z a zero 
bias conductance peak emerges (curves b and c). The resulting ZBCP is broad in 
contrast to the spin-singlet af-wave case due to the fact that the position of the ABS 




Fig. 3 (Color online) (a): Tunneling conductance for a chiral p-wave superconductor, a: Z = 0, b: 
Z = 1 and c: Z = 5. (b): Tunneling conductance for a chiral p-wave superconductor in the presence 
of a magnetic field for Z = 0. a: H = 0, b: H = Q.2H and c: H = QAH , d: H = -0.2H and e: 
H = -QAH Q . 
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depends on the injection angle according to Eq. (TBI [8 1. The presence of the ABS 
has been confirmed by tunneling experiments I1121I13II . 

Next we consider the situation where a magnetic field H is applied perpendic- 
ular to the two-dimensional plane, which induces a shielding current along the in- 
terface. When the penetration depth for the chiral p-wave superconducting mate- 
rial is much longer than the coherence length, the vector potential can be approx- 
imated as A(r) = (0,A v (x),0) with A y (x) — —A m Hexp(—x/A. m ), where Am is the 
penetration depth. In the following we consider the situation where Landau level 
quantization can be neglected. Then the quasiclassical approximation can be used. 
The applied magnetic field shifts the quasiparticle energy E to E + HAo&m(p /Hq 
with Hq = h/ (2e7T 2 ^A, n ) and | = h 2 ^ / (nmAo) [29|. The resulting tunneling con- 
ductance for various magnetic fields is plotted in Fig. [3jb). As is seen, Oj(E) is 
enhanced for positive H, while it is reduced for negative H. This can be roughly 
understood by looking at the bound state levels. In the presence of H, the bound 
state energy can be expressed by 

E b (6)=A (l -H/Hq) sin ~4o(l -H/H Q )k y /k P . (15) 

The contribution of the Andreev bound state to the conductance enters via a term 
8(E —Eb(0)), which is proportional to 1 /|d£'b(0)/d0|. It is clear that the slope of 
the dispersion around = is reduced for positive H, leading to an enhancement 
of the numerator in Eq.[lO] around = 0, where the bound states are close to zero 
enery. On the other hand, for negative H, the height of the ZBCP is reduced since 
the slope of the curve of Et> around zero energy becomes steeper [ 10] . 

In p-wave superconductors self-consistency of the order parameter and impurity 
effects can be of importance. In Fig.|4]we show self-consistent order parameters and 
Andreev spectra at a surface of a layered p-wave superconductor. In addition to the 
bulk k x + iky component a subdominant k x — ik y component is stabilized within a 
few coherence lengths = vp/2nkBT c ) near the surface. The full lines are results 
assuming a mean free path of I = 10<^o everywhere. When replacing the mean free 



Fig. 4 (Color online) Self- 
consistent order parameter 
near a surface of a layered 
p-wave superconductor. Full 
lines are for mean free path 
^mfp = lOiijo everywhere; 
dashed lines are for a short- 
ened ^mfp = 03^o in the gray 
shaded region. The calcu- 
lations are for T = 0.ir c . 
Insets: point contact spectra 
for fully transparent interface 
(bottom) and tunneling con- 
ductace (O = 0.05) (top). 
After Ref. j30). 
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path in a surface layer (gray shaded region in Fig. |4j) by £ = 0.3<i;o, we obtain the 
results shown as dashed lines. In contrast to the first case, for the second case both 
order parameter components are strongly suppressed near the surface. The presence 
of an increased scattering in a surface layer also modifies the form of point contact 
spectra and the tunneling conductance as seen in the insets of Fig. |4] In constrast 
to the surface density of states which for a clean surface is constant in energy, the 
tunneling conductance shows a broad peak similar as in Fig. [3] which is however 
reduced in height for a self consistent order parameter J3T|. 

Finally, we would like to comment that the above edge state is topologically 
equivalent to that of a quantum Hall system. In a quantum Hall system it is estab- 
lished that the edge channel supports the accurate quantization of the Hall conduc- 
tance Oh, which is related to a topological integer [32, 33]. In the edge state of a 
chiral /?-wave superconductor, such a topological number can be also defined l34ll . 
For this case, the edge state is topologically protected by the bulk energy gap Aq. 
The topological properties of the electronic states have been attracting intensive in- 
terest in condensed matter physics. In section l2.3.3l we will return to this question 
in connection with non-centrosymmetric superconductors. Before that, we discuss 
in the following section theoretical predictions for the Andreev conductance spectra 
for non-centrosymmetric superconductors. 



2.3 Andreev conductance in non-centrosymmetric superconductors 

Non-centrosymmetric superconductors such as CePt3Si are a central topic of current 
research [35,36]. Two-dimensional non-centrosymmetric superconductors are ex- 
pected e.g. at interfaces and/or surfaces due to a strong potential gradient. An inter- 
esting example is superconductivity at a LaA103/SrTi03 interface |37, 38 1. In non- 
centrosymmetric materials spin-orbit interaction becomes very important. Frigeri et 
al. IT361 have shown that the (p x ± \p y ) -pairing state has the highest T c within the 
triplet-channel in CePt^Si. It has been shown that singlet (s-wave) and triplet (p- 
wave) pairing is mixed, and several novel properties related to that mixing, such as 
a large upper critical field beyond the Pauli limit, have been focused on [36|. On 
the other hand, a pure (p x ± ip v )-pairing state has been studied as a superconduct- 
ing analogue of a quantum spin Hall system [39|. Therefore, it is an important and 
urgent issue to study the spin transport properties of the NCS superconductors from 
a topological viewpoint. 

In this section, we discuss charge and spin transport in non-centrosymmetric su- 
perconductors |40 1 . We concentrate on non-centrosymmetric superconductors with 
time-reversal symmetry, where a spin-triplet (p x ± i/?j,)-wave and a spin-singlet s- 
wave pair potential can mix with each other, similar as discussed in the last section. 
We show that when the amplitude of the (p x ±ip y )-wave component is larger than 
that of the s-wave component, then the superconducting state belongs to a topolog- 
ically nontrivial class analogous to a quantum spin Hall system, and the resulting 
helical edge modes are spin current carrying Andreev bound states that are topo- 
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logically protected. Below, we study Andreev reflection PTll at low energy, which 
is determined mostly by the helical edge modes, and find the spin polarized current 
flowing through an interface as a function of incident angle. When a magnetic field 
is applied, even the angle-integrated current is spin polarized. 

2.3.1 Andreev bound states 

We start with the Hamiltonian of a non-centrosymmetric superconductor 

H(k) A(k) \ 
-A*(-k) -H*(-k) J 

with H{k) = | t + g(k) a, g(k) = X (xk y - yk x ), % k = Ti 2 k 2 /(2m) - p. Here, jj,, 
m, <7 and A denote chemical potential, effective mass, Pauli matrices and coupling 
constant of Rashba spin-orbit interaction, respectively [ 36 ] . The pair potential A (k) 
is given by 

A(k) = [d(k)-a + Y(k)]ia y . (16) 

We choose (p x ±ip y ) with d(k) = A p (xk y —yk x )/ \ k | for the spin-triplet component 
ll36ll . and y(k) = A s with A p > and A s > 0. The superconducting gaps A\ = A p + 
A s and A2 =\A p — A s \ open for the two spin-split energy bands, respectively, in 
the homogeneous state [42|. As we will show below, surface states are crucially 
influenced by the relative magnitude between A p and A s . 

Let us consider the wave functions, focusing on those for ABS localized at the surface. Consider 
a two-dimensional semi-infinite superconductor forx > where the surface is located atx = 0. The 
corresponding wave function is given by 1 43 1 

%(x) =e k ? y [ci^ 1 e i *i^ + c 2 V'2e _i? ^ + C3y3e i «2^ + C4^4e^ ,( l , 

4 = *£ ± (kj/l&JiEi-AjyiV+Vtfi/m), (17) 

with j = {1,2}, and kf x = kj x = kj x for | k y |< kj and kj x = — kj x = kj x for | k y > kj. Here, k\ and 
ki with 

k m = zprnA /h 2 + sJ(mX/h 2 ) 2 +2mn/h 2 (18) 
are the Fermi momenta of the small and large Fermi surface, respectively (the upper sign holds for 
k\ ), and kj x denotes the x component of the Fermi momentum kj, with kj x = \Jk 2 — k 2 . The wave 
functions are given by 



Yl = 



1 «1 
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( V ' 
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/ 112 \ 
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— iaf 'w 
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with y= sgn(4 p — A s ). In the above, 

Uj = ME + ^E^A 2 )/2E, VJ = J(E-^E^A])/2E. (19) 
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Here we have introduced ai = (fet — \ky)/k\ , a-i = (k^ — ik y )/k2, &\ = (— — iky)/ki, and 
a,2 = ( — k^ x — iky ) /&2 ■ E is the quasiparticle energy measured from the Fermi energy. By postulating 
ft; (x) = at x = 0, we can determine the ABS. 

The bound state condition can be expressed by 

(4f-£2)(4|-£2) = I-| (£2 + 74,42), (20) 
f^l±M for | 02 |< 0C 

£= <^ C S 2[^(ei-»2)] 1 2 — (21) 

U for C <| 2 |< tc/2, 

with £ < 1, cos 0i = k\ x /k\ and cos 02 = kixjki. The critical angle C is defined 
as arcsin(&i/&2)- For A = 0, Eq. ( 120b reproduces the previous results 11421 . As seen 
from Eq. d20l ). a zero energy ABS is only possible for | 02 |< C and 7=1, i.e. 
4 P > -4 S . This ABS corresponds to a state in which a localized quasiparticle can 
move along the edge. The energy level of this edge state depends crucially on the 
direction of the motion of the quasiparticle. The inner gap edge modes are absent for 
large magnitude of k y , i.e. 02- In this case, k\ x becomes a purely imaginary number 
due to the conservation of the Fermi momentum component parallel to the surface. 
The parameter regime where the edge modes survive is reduced with increasing X. 
However, as far as we concentrate on normal injection, the edge modes survive as 
midgap ABS (2[5| irrespective of the strength of A. If we focus on the low energy 
limit, the ABS energy can be written as 

/ A}\k x +k 2 

E = ±A "{ l -^ p )hi k - (22) 

with A s < A p for any X with small magnitude of k y . For A s > A p , the ABS vanishes 
since the value of right hand side of Eq. ( f20b becomes negative, due to the negative 
sign of 7 for | E |< A\ and \E \< Aj_. It should be remarked that the ABS under 
consideration does not break time reversal symmetry, since the edge currents carried 
by the two partners of the Kramers doublet flow in opposite directions. Thus they 
can be regarded as helical edge modes, with the two modes related to each other by 
a time reversal operation. 



2.3.2 Charge and spin conductance 

Now we turn to transport properties governed by the ABS in NCS superconduc- 
tors [44-46|. First, we point out that the spin Hall effect, i.e., the appearance of 
the spin Hall voltage perpendicular to the superconducting current, is suppressed 
by the compressive nature of the superconducting state by the factor of (^F^m) 2 
(kp: Fermi momentum, Am: penetration depth) [34|. Instead, we will show below 
that spin transport through the junction between a ballistic normal metal at x < 
and a NCS superconductor, i.e., through a N/NSC junction, can be enhanced by 
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the Doppler effect during Andreev reflection. The Hamiltonian Mfa of N is given 
by putting A(k) = and A = in We assume an insulating barrier at x = 0, 
expressed by a delta-function potential U8(x). 

The quantities of interest are the angle resolved spin conductance fs(0) and 
charge conductance fc(0) defined by ll47ll 

/s( ) = \ [ E *p(I a °-P I' - I Vp I')] cos > f 23 ) 

a.p 

fc(0)= [l + iE(l"^pl 2 -|^pl 2 )] cos0 ' C24) 

L a.p 

where s p = +(— )1 for p =t (4), and 9 denotes the injection angle measured 
from the normal to the interface. Here, b ap and a a , p with a,p € {t>4} are spin- 
dependent reflection and Andreev reflection coefficients, respectively. 

These coefficients are determined as follows. The wave function for spin o~ in the normal metal 
tffo (x) is given by 

f^(jr)=exp(i^ Fv y)[(l^ CT + £ a a . p y ap )exp(ik R x) + £ fo (J .pVt p exp(-i/c Fl x)] (25) 
P=U P=U 

with T Yn = T Ybt = {10,0,0), T Vil = T Vb , = (0, 1 , 0, 0), T Vat = (0, 0, 1 , 0), and T Vai = (0, 0, 0, 1 ) . 
The corresponding 1%(x) is given by Eq. (T7J. The coefficients a a p and & CT ,p are determined by 
postulating the boundary condition ffr(0) = f s (0), and M Sx W s (0)'- hv Nx ^(0) = -2;f/f 3 f / s (0) 
with fivs(N)i = dHs(N)/dkx' an< ^ me diagonal matrix T3 given by % = diag(l, 1,-1,-1). 

The resulting angle averaged charge conductance (tunneling conductance) is 
given by 

at = cjt = ( /^/c(e)<ifl) / ( f_^ 2 hc{B)dd^ . (26) 

We plot in Fig.[5]the charge conductance by changing the ratio of 4 S / in the pres- 
ence of the splitting of the Fermi surface B5I . For 4 4 < Ap, OT(eV) has a ZBCP due 
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to the presence of the helical edge modes (curves a and b in Fig. |5j. For A s = A p , 
due to the closing of the bulk energy gap, the resulting Gj(eV) is almost constant. 
For A s > A p , Oi(eV) has a gap like structure similar to spin-singlet s-wave super- 
conductor. 

Next, we focus on the spin conductance. First we consider a pure (p x ±ip y )- 
wave state. In Fig. [6] the angle resolved spin conductance is plotted as a function of 
injection angle 9 and bias voltage V with E = eV . Note here that the k y is related 
to 9 as k y = kpsinQ. It is remarkable that the spin conductance has a non zero 
value although the NCS superconductor does not break time reversal symmetry. The 
quantity fs(9) has a peak when the angle 9 or k y corresponds to the Andreev bound 
state energy E in the energy dispersion. With this condition, the spin-dependent 
Andreev reflection occurs to result in a spin current. Besides this property, we can 
show that fs{9) = — /s(— 0) is satisfied. By changing the sign of eV, fs{9) changes 
sign as seen in Fig. |6ja). Next, we look at the case where an .s-wave component 
coexists. We calculate the spin conductance similar to that for the pure (p x ± ip y )- 
wave case. For A s < A p , where helical edge modes exist, fs(9) shows a sharp peak 
and/s(0) = — /s(— 9) is satisfied [see Fig.|6jb)]. These features are similar to those 
of the pure (p x ± ip v )-wave case. On the other hand, for A s > A p , where the helical 
edge modes are absent, sharp peaks of fs{9) as shown in Fig.|6]are absent. 

We have checked that there is negligible quantitative change, i.e., less than 0.5% 
change of the peak height, by taking the A = limit compared to Fig. [6] In this limit, 
for the pure (p x ± ip y )-wave state, fs(9) is given simply as follows 

-8/?D 2 sin20sin2<pcose 
| 4(sin 2 9 - sin 2 <p) +D[2cos20 - (1 +fl)exp(-2i<p)] | 2 




Fig. 6 (Color online) Angle resolved spin conductance for Z = 5. a: eV = 0. \A p , b: eV = — 0. \A p 
and c: eV = 0.64 p with Xkf = 0.1^. (a) pure (p x ±ipj,)-wave case with A s = 0; (b) A s = Q3A P . 
[from Fig. 2 Phys. Rev. B 79, 060505(R) (2009).] 
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for E < A p and fs(9) — for | E |> A p with sin<p = E/A p The transparency 
of the interface D is given as before by 4cos 2 0/(4cos 2 9 +Z 2 ), with the dimen- 
sionless constant Z = 2mU /h 2 kp. The magnitude of fs(9) is largely enhanced at 
E = zkApSinO corresponding to the energy dispersion of the ABS. The origin of 
the nonzero fs(0) even for X — is due to spin-dependent Andreev bound states. 
We have checked that even if we take into account the spatial dependence of the 
(p x ± if> v )- wave P a ir potential explicitly, the resulting /s(0) does not qualitatively 
change fi4l . 

Summarizing these features, we can conclude that the presence of the helical 
edge modes in NCS superconductors is the origin of the large angle resolved spin 
current through normal-metal/NCS superconductor junctions. However, the angle 
averaged normalized spin conductance becomes zero since fs{9) = —fs(—9) is 
satisfied. 

Magnetic field offers an opportunity to observe the spin current in a more acces- 
sible way, where the time reversal (T) symmetry is broken by the shielding current 
at the interface. Here we consider the angle averaged normalized spin conductance 
Cs and charge conductance afc as a function of magnetic field. The spin conductance 
is given by ||47l 

aS = ( /l// 8 ^^ ) / ( lTi2 hc{d)dd ) ' ^ 

where /nc(^) denotes the angle resolved charge conductance in the normal state 
with A p = A s = 0. We consider a magnetic field H applied perpendicular to the two- 
dimensional plane, which induces a shielding current along the normal-metal/NCS 
superconductor interface. When the penetration depth of the NCS superconductor 
is much longer than the coherence length, the vector potential can be approximated 
as described in section 12.21 As in the case of a chiral p-wave superconductor, the 
applied magnetic field shifts the quasiparticle energy E to E + HA p smO /Hq. For 
typical values of | ~ 10 nm, Am ~ 100 nm, the magnitude of Hq is of the order of 
0.2 Tesla. The order of magnitude of the Doppler shift is given by HA p /Hq. Since 
the Zeeman energy is given by jX^H, the energy shift due to the Doppler effect is 
by a factor kpX m larger than that due to the Zeeman effect. Thus, we can neglect 
the Zeeman effect in the present analysis. This is in sharp contrast to quantum spin 
Hall systems where the Zeeman effect is the main effect of a magnetic field, which 
opens a gap in the helical edge modes and modulates the transport properties [48 1 . 
The enhanced spin current due to Doppler shifts is specific to the superconducting 
state, and is not realized in quanum spin Hall systems. 

2.3.3 Topological aspects 

We now focus on the topological aspect of non-centrosymmetric superconductors. 
Recently, the concept of the quantum Hall system has been generalized to time- 
reversal (T) symmetric systems, i.e., quantum spin Hall systems Il49| - |5T1 . A quan- 
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turn spin Hall system could be regarded as two copies of a quantum Hall system, 
for up and down spins, that are characterized by opposite chiralities. In the generic 
case, however, a mixture of up and down spins occurs due to spin-orbit interac- 
tion, which necessitates a new topological number to characterize a quantum spin 
Hall system 149115 II . In quantum spin Hall systems, there exist helical edge modes, 
i.e., time-reversed partners of right- and left-going one-dimensional modes. This 
has been experimentally demonstrated for the quantum well of the HgTe system by 
measurements of the charge conductance J48). 

As shown in Fig. [6] to discuss the topological nature of the helical edge modes, 
it is sufficient to consider the pure (p x ± i/?,,)-wave state. Here, we give an argument 
from the viewpoint of the Zi (topological) class ||49l , why the superconducting state 
with A p > A s has an Andreev bound state. We commence with a pure (p x ± ipy)- 
wave state without spin-orbit interaction, i.e. A = 0. The spin Chern number llBTI for 
the corresponding Bogoliubov-de Gennes Hamiltonian is 2. Turning on A adiabati- 
cally leaves the time reversal T-symmetry intact and keeps the gap open. Upon this 
adiabatic change of A, the number of the helical edge mode pairs does not change. 
The reason is that this number is a topological number and consequently can only 
change by integer values. We now increase the magnitude of A s from zero. As far as 
A p > A s is satisfied, the number of helical edge modes does not change. However, if 
A s exceeds A p , the helical mode disappears. In this regime, the topological nature of 
the superconducting state belongs to a pure s-wave state with A = 0. It is remarkable 
that just at A s = A p one of the two energy gaps for quasiparticles in the bulk closes. 
At precisely this point a quantum phase transition occurs. 

In the following, we discuss the pure (p x ±ip y )-wa\e case in more detail. In Fig. 
[7] the spin conductance Os an d charge conductance Gq normalized by the charge 
conductance in the normal state are plotted. It should be noted that Os becomes 




Fig. 7 (Color online) Angle averaged spin conductance and charge conductance as a function of eV 
with bias voltage V with Xkp = O.ljU. a: H = 0, b: H = -0.2H , c: H = 0.2H , and d: H = -0AH . 
Curves b and c of the right panel are identical, [from Fig. 2 Phys. Rev. B 79, 060505(R) (2009).] 
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nonzero in the presence of a magnetic field H (see curves b, c and d), since f$ (9) is 
no more an odd function of 9 due to the imbalance of the helical edge modes. For 
A = 0, the corresponding helical edge modes are given by E = A p (1 —H /Ho) sin 9 
and E = — A p (l + H/Hs) sin 0. As seen from the curves b and c, the sign of Os is 
reversed when changing the direction of the applied magnetic field. On the other 
hand, the corresponding charge conductance has different features. For H = 0, the 
resulting line shape of afc is the same as that for a chiral p-wave superconductor 
(see curve a of right panel) [42 43 45 1. As seen from curves b and c in the right 
panel, afc does not change with the direction of the magnetic field H. 

In summary, we have clarified the charge and spin transport properties of non- 
centrosymmetric superconductors from the viewpoint of topology and Andreev 
bound state. We have found an incident angle dependent spin polarized current 
flowing through the interface. When a weak magnetic field is applied, even the 
angle-integrated current is largely spin polarized. In analogy to quantum spin Hall 
systems, the Andreev bound states in non-centrosymmetric superconductors corre- 
sponds to helical edge modes. Andreev reflection via helical edge modes produces 
the enhanced spin current specific to non-centrosymmetric superconductors. 



3 Quasiclassical Theory of Superconductivity for 
Non- Centrosymmetric Superconductors 

3.1 Quasiparticle Propagator 

Electronic quasiparticles in normal Landau Fermi liquids are restricted in phase 
space to a region that comprises only a small part of the entire electronic phase 
space 15211531 . It consists of a narrow (compared to the Fermi momentum p F ) shell 
around the Fermi surface, and a small (compared to the Fermi energy E-p) region 
around the chemical potential. Quasiparticles are characterized by their spin and 
charge, and their group velocity is the Fermi velocity, Vp(p F ). Quasiclassical theory 
is the appropriate framework to describe such a system. It consists of a system- 
atic classification of all interaction processes according to their relevance, i.e. their 
smallness with respect to an expansion parameter SMALL I54H581 . This expansion 
parameter assumes the existence of a well defined scale separation between a low- 
energy scale and a high-energy scale. 

Superconducting phenomena are governed by the low-energy scale. That means 
that the energy scales determined by the energy gap A and the transition tempera- 
ture T c are small. In contrast the energy scales determined by the Fermi energy Ep 
or the Coulomb repulsion Uc are large energies. Disorder can be described within 
quasiclassical approximation as long as the energy associated with the scattering 
rate, h/x, is classified as a small energy. A systematic classification shows that a 
consistent treatment of disorder requires the f -matrix approximation. Localization 
effects due to disorder are beyond the leading order precision of quasiclassical the- 
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ory. Associated with the energy scales are small and large length scales. For example 
the superconducting coherence length = hvp/27lks,T c , and the elastic mean free 
path I = VfT are large compared to the lattice constant a and the Fermi wave length 
A F = h/p F . 

This separation in energy and length scales is associated with a low-energy region 
in phase space, that includes low quasiparticle energies ~ A,k^T, and a momentum 
shell around the quasiparticle Fermi momentum p F of extend 8p ~ A / 1 Vf (p F ) \ ■ The 
phase space volume of this low-energy region, divided by the entire phase space vol- 
ume, is employed for a systematic diagrammatic expansion of a Dyson series within 
a path-ordered Green's function technique (e.g. Matsubara technique for the Mat- 
subara path, Keldysh-Nambu-Gor'kov technique for the Schwinger-Keldysh path). 
Within the framework of Green's function technique, all diagrams in a Feynman 
diagrammatic expansion can be classified according to their order in this expansion 
parameter, which is denoted as SMALL. The leading order theory in this expan- 
sion parameter is called the "Quasiclassical Theory of Metals and Superconduc- 
tors" munio). 

The possibility to define a quasiparticle Fermi surface around which all quasi- 
particle excitations reside is a requirement for the quasiclassical theory to work. Its 
presence ensures that the Pauli principle is still effective in placing stringent kinetic 
restrictions on the possible scattering events. It is essential to note that such a defini- 
tion need not be sharp, i.e. the theory is not restricted to normal Fermi liquids with 
a jump in the momentum distribution at zero temperature. Thus, the theory includes 
superconducting phenomena as well as strong coupling metals. It is convenient to 
introduce a local coordinate system at each momentum point of the Fermi surface 
p F , with a variation along the surface normal, i.e. in direction of the Fermi veloc- 
ity Vp(pp), that is determined by a variable £, p (this variable is zero at the Fermi 
momentum), and a tangential variation along the Fermi surface at constant £, p . A 
consistent approximation requires to consider the Fermi velocity constant across 
the low-energy momentum shell, and thus the local coordinate system stays an or- 
thogonal system as long as £, p varies within this momentum shell, and furthermore, 
^ p stays small within this momentum shell. The coordinate £, p around each Fermi 
surface point p F varies then approximately as ^ p m v F (p — p F ). 

The quasiclassical theory is obtained by defining quasiparticle propagators for 
the low-energy regions of the phase space, and in combining all diagrams involv- 
ing Green's functions with their variables residing in the high-energy regions into 
new effective high-energy interaction vertices. This process of integrating out high- 
energy degrees of freedom is highly non-trivial and must be solved by microscopic 
theories. In the spirit of Fermi liquid theory it is, however, possible to regard all 
high-energy interaction vertices as phenomenological parameters of the theory. In 
quasiclassical approximation they do not depend on any low-energy variables as 
temperature or superconducting gap, and they do not vary as function of ^ p as long 
as <^ p stays within the momentum shell that harbors the quasiparticle excitations. 
However they do depend in general on the position of the Fermi momentum on the 
Fermi surface. 



Surface properties of non-centrosymmetric superconductors 



17 



In addition to introducing new effective interaction vertices the above procedure 
also introduces a quasiparticle renormalization factor a 2 (p) ~ 1/Z(p), that is due 
to the self energies of the low-energy quasiparticles moving in the background of 
the high energy electrons. This renormalization leads to a modification of the quasi- 
particle Fermi velocity compared to the bare Fermi velocity of the system, and to a 
deformation of the quasiparticle Fermi surface compared to the bare Fermi surface. 
It also determines quasiparticle weight as the residua of the quasiparticle poles in 
the complex energy plane. 

One has to keep these remarks in consideration when including additional inter- 
action, like spin-orbit interaction or exchange interaction, in a quasiclassical theory. 
First, it is important to decide if this interaction is going to be treated among the low- 
energy terms or among the high-energy terms. Depending on that one obtains two 
different quasiclassical theories, that cannot in general adiabatically be connected 
with each other. Going from one limit to the other includes the un-dressing of all ef- 
fective interaction vertices and of the quasiparticles, and re-dressing with new types 
of effective interaction vertices and self energies. Importantly, this dressing leads to 
strongly spin dependent effective interactions and quasiparticle renormalizations in 
one limit, and to leading order spin-symmetric interactions and quasiparticle renor- 
malizations in the other limit. The former case, when spin-dependent interactions 
are included in the high-energy scale, leads to a complete reorganization of the 
Fermi surface geometry, with in general new spin-dependent quasiparticle energy 
bands. In this case, it is not sensible anymore to keep the spin as a good quantum 
number, but it is necessary to deal directly with the representation that diagonalizes 
the energy bands including the spin-dependent interaction. In the case of a strong 
exchange energy this leads to exchange split energy bands, and in the case of strong 
spin-orbit interaction this leads to helicity bands. 

The basic quantities in the theory are the quasiparticle Fermi surface, the quasi- 
particle velocity, and quasiparticle interactions. Here we give a short sketch of how 
they enter the theory. The bare propagator (without inclusion of exchange interac- 
tion or spin-orbit coupling) in the quasiparticle region of the phase space has the 
general structure 



G (°) i n e ) - ^ 



(28) 



where £ ^°\p) is the bare energy dispersion of the energy band (measured from the 
electrochemical potential of the electrons). It does not include electron-electron in- 
teraction effects yet, and thus determines a bare Fermi surface that does not coincide 
with the quasiparticle Fermi surface defined below. The quantum number a labels 
the spin. The leading order self energy is solely due to coupling of low-energy elec- 
trons (superscript L) to high-energy electrons (superscript H), and consequently the 
corresponding self energy, Z^ H \ must be classified as a pure high-energy quantity. 
In general, when either exchange interaction or spin-orbit coupling are large energy 
scales, this self energy contribution will be spin-dependent (and will ultimately lead 
to new, spin-split energy bands as explained below). The self energy is, how- 
ever, slowly varying in energy on the low-energy scale and thus can be expanded 
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around the chemical potential, 

^S ) (p^)=^S ) ^°) +e ^S ) ^ e )L = o + ^ e2 )- (29) 

The second term can be combined with energy e into a renormalization function 



e-ed e Z§j(p,0)=Z < §(p)e, (30) 

such that to leading order the Dyson equation for the low-energy propagator G^) 
reads 

{zg(p)e [^\p)d a p+Z^(j,,0}-L^{p,e)}^G^{p,e) = 8 ar , (31) 

where includes all self energy terms of order SMALL . Here, and in the follow- 
ing, summation over repeated indices is implied. The (g) sign accounts for possible 
spatial or temporal inhomogeneities, in which case it has the form of a convolution 
product in Wigner representation (see Ref. [54] for details). The equation holds in 
this form either in Matsubara or in Keldysh representation (in which case all quanti- 
ties are 2x2 matrices in Keldysh space [67]). Low-energy excitations reside in mo- 
mentum regions differing considerably from that for the bare propagators Eq. (1281 1. 

The quantities Z K J(j>) and E a p (p,0) can be defined such that they have real 
eigenvalues. The next step is to eliminate the high-energy renormalization factor 
Z a p(p) from the low-energy theory. This is done with the help of quasiparticle 
weight factors a a p(p), that are the solution of 

aay(p)Zyy [p)ayp (p) = 8 a p . (32) 

They exist as long as Z^y has non-zero eigenvalues. Then we can define the quasi- 
particle Green's function G^p as the solution of 

aar(p)G^ P) ayp(p) = G§>(p,e), (33) 

which exists under the condition that a a y has non-zero eigenvalues (i.e. the quasipar- 
ticle weights are non-zero; otherwise the quasiparticle approximation breaks down). 
It fulfills the Dyson equation 

[e -$Q p (p)-ZQ p (p,e)] ap ®Gf Y P \p,e) = 8 ay (34) 
with the quasiparticle dispersion 

$af \P) = a *r(P) (^ 0) (p)8 r y+Eff(p,0)) ayp(p) (35) 
and the quasiparticle self energies 
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E aB (P^)= a ay(p)KJ(P,£)afp(p). 



(36) 



The effective (renormalized by high-energy processes) interactions vertices for the 
low-energy propagators, G^, which enter the diagrammatic expressions for the 
quasiparticle self energy, have the general structure Vp l ^p n (£i,Pi;...;£ n ,p n ). In 
leading order the energy dependence of these vertices can be neglected near the 
chemical potential, i.e. the arguments can be restricted to the chemical potential. 
Furthermore, instead of working with G^j and V^...^ the common and com- 
pletely equivalent description in terms of the above defined quasiparticle propa- 
gators, G^p \ and renormalized quasiparticle interactions, vjf^ P ^, given by 



y ft (Pl —Pn) = fl A«(Pi) • ••«/*.« (Pn) v m...l}> (0,Pi; • • • ;0,p„) (37) 



can be used. 

It is important to note that the quasiparticle self energies can be written down as 
functionals of the quasiparticle Green's functions only in leading order in the expan- 
sion in SMALL , which is the order relevant for the quasiclassical approximation. In 
this case, the quasiparticle weights have disappeared from the theory and cannot in 
principle be determined from low-energy processes that only involve quasiparticle 
dynamics. They must be obtained from a microscopic theory by considering high 
energy scattering processes, which is beyond the quasiclassical approximation. 

It is obvious, that the appearance of the quasiparticle renormalization factors 
renders all self energies and interactions non-diagonal in spin unless spin-dependent 
interactions are small enough the be omitted from the high-energy quantities. From 
the above expressions one obtains the quasiparticle Fermi surfaces by diagonalizing 
the quasiparticle dispersion 




(38) 



with band index A, and solving the equation 



& qp) (p) = o^p = p x f . 



(39) 



The corresponding quasiparticle Fermi velocity is then given by 




(40) 



In the band diagonal frame, the quasiparticle propagator is given by, 



{[e-^ p (p)]8 XXl (/>,£)} ®Gg?(p,e) = 8 XX , 



(41) 



where the self energy (and all interactions in the self energy expressions) must be 
transformed accordingly, e.g. 
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E^{p,e) = U p , a E^{p,e)u; pi , (42) 

In the next section this procedure is carried out for the case of a strong spin-orbit 
interaction, e.g. appropriate for some non-centrosymmetric materials. 



3.2 Spin-orbit interaction and Helicity representation 

As discussed in the introductory chapter of this book, for treating a non-centrosym- 
metric material it is convenient to perform a canonical transformation from a spin 
basis with fermion annihilation operators a ka for spin a =t,4 to the so-called helic- 
ity basis with fermion annihilation operators c k x for helicity X = ±. This canonical 
transformation diagonalizes the kinetic part of the Hamiltonian, 

■^ = 1 E [S(t)+»(*)-o)a/i]oU = II^U- («) 

k ap=f,l k l=± 

Here, % (Jfc) is the band dispersion relative to the chemical potential in the absence of 
spin-orbit interaction, g(k) is the spin-orbit pseudovector, which is odd in momen- 
tum, g(—k) = —g(k), and a is the vector of Pauli matrices. The resulting helicity 
band dispersion is 

Uk) = m±\g(k)\- (44) 

As is easily seen, spin-orbit interaction locks the orientation of the quasiparticle spin 
with respect to its momentum in each helicity band. The Hamiltonian, Eq. d43l l. is 
time reversal invariant, however lifts the spin degeneracy. 

It is convenient to introduce polar and azimuthal angles for the vector g, defined 
by {g x ,gy,8z} = l«|{sin(0g)cos((j!)g),sin(0g)sin((j!> g ),cos(0 g )} (whereO < 9 g < n). 
In terms of those, the transformation from spin to helicity basis, U k i a , is defined 
by 061 

/ cos(0 g /2) sin(0j/2)e-''M _ v 
UkXa- ^_ s in(e j /2)e , 'ft COs(0 g /2) J 1 c kX ~ ^kXa^ka- (45) 

Obviously, Zap U kXa [g (k) ■ O a 0J^ = \g(k)\a$,. 

For the superconducting state the Nambu-Gor'kov formalism is appropriate [66]. 
The Nambu spinor, A k = (at[,o k i,c^_ k *.,a_ k ^) T transforms under the above canon- 
ical transformation into the helical object C k = (cjt+,Cjt-,c^ t+ ,c^ t ) r , where 

C k = U k A k , #*=(o*i/°*J- (46) 
Correspondingly, one can construct 4x4 retarded Green's functions in spin basis, 
G { k!k 2 (tht2) = -ie(t 1 -t 2 )({A kl (t 1 ),Al 1 (t2)})je, (47) 
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and in helicity basis, 

<W'i,f2) = -t2)({C kl (h),cl 2 (t 2 )})je,= t) kl G^ k2 (ti,t2)Ul 2 , (48) 

where A(t) and C(t ) are Heisenberg operators, the braces denote an anticommuta- 
tor, (■ ■ - is a grand canonical average, and 9 is the usual Heaviside step func- 
tion. Analogously, advanced, Keldysh, and Matsubara propagators can be defined 
in helicity representation. For dealing with superconducting phenomena it is often 
convenient to introduce Wigner coordinates, 

G(k,R,e,t) = l(dq)(dT)^ R+ ^G k+u _,_(t + ^t-^). (49) 

From here, one can proceed along different lines. Either, the Dyson equation for the 
full Gor'kov Green's functions is solved, which is equivalent to the Bogoliubov-de 
Gennes description in wave function techniques. Or, the quasiclassical approxima- 
tion is employed, that is equivalent to the Andreev approximation in wave function 
language. In the following section we will adopt the second line. 



3.3 Quasiclassical Propagator 

In the following, the quasiclassical theory of superconductivity [54,59 -65] will be 
employed to calculate electronic transport properties across interfaces with non- 
centrosymmetric superconductors. This method is based on the observation that, 
in most situations, the superconducting state varies on the length scale of the super- 
conducting coherence length §q = hvp/27zkgT c . The appropriate many -body Green's 
function for describing the superconducting state has been introduced by Gor'kov 
[66], and the Gor'kov Green's function can then be decomposed in a fast oscillat- 
ing component, varying on the scale of 1 / kp, and an envelop function varying on 
the scale of <^q. The quasiclassical approximation consists of integrating out the fast 
oscillating component for each quasiparticle band separately: 

g(p%,R,e,t) = J d# f 3 G (QP Hp,R,e,t) (50) 

where a "check" denotes a matrix in Keldysh-Nambu-Gor'kov space, l67l a "hat" 
denotes a matrix in Nambu-Gor'kov particle-hole space, = Vp (p — Pp), and T3 
is the third Pauli matrix in particle-hole space. 

The quasiclassical Green's function obeys the transport equation [59 60 1 

ihv F -V R g + [ef 3 -A-h,g} = b. (51) 

Here, e is the quasiparticle energy, A is the superconducting order parameter and h 
contains all other self-energies and external perturbations, related to external fields, 
impurities etc. The notation o combines a time convolution with matrix multiplica- 
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tion, and [•,•]<, denotes the commutator with respect to the o-product. Equation d5TT l 
must be supplemented by a normalization condition that must be obtained from an 
explicite calculations in the normal state 05911681 , 

gog=-W. (52) 

From the knowledge of g one can calculate measurable quantities, e.g. the current 
density is related to the Keldysh component of the Green's function via 

/ft F 
— Tr<v F f 3 | K CpM,e,f)>, (53) 

where q = —\e\ is the electron charge, and (• • •) denotes a Fermi surface average, 
which is defined by 

and Tr denotes a trace over the Nambu-Gor'kov matrix. 



3.3.1 Case of weak spin-orbit splitting 

In the case of weak spin-orbit splitting the quasiclassical propagator can be obtained 
in either spin or helicity representation. It is possible then to define a common Fermi 
surface p F for both spin bands or, equivalently, both helicity bands. This case applies 
when \g(p-p) | <C Ep for any Fermi momentum p F , where Ep is the Fermi energy (in 
addition to the condition that the superconducting energy scales {ksT c and the gap A 
are much smaller than Ep). Under these circumstances quasiparticles with different 
helicity but with the same k = k/\k\ propagate coherently along a common clas- 
sical trajectory over distances much longer than the Fermi wavelength. The trans- 
port equation is the usual Eilenberger equation modified by a spin-orbit interaction 
term 154ll69l 

3tVF-Vj{|+[ef3-2-v S o,|] = 6 (55) 

with normalization g og — — 7T 2 1. Here, in helicity basis vso — \&k v \ O^K and in spin 
basis vso — gk F ' **^3' w i tn 

d= (oaO = (o-a(W))- (56) 

The velocity renormalization of order \g\/Ep <C 1 can safely be neglected. The 
quasiparticle trajectories are doubly degenerate in either spin or helicity space, and 
coherent mixing between spin states or between helicity states can take place. 
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3.3.2 Case of strong spin-orbit splitting 

In the case of strong spin-orbit splitting the only possible representation for quasi- 
classical theory is the helicity representation. In this case, the spin-orbit interaction 
does not appear anymore as a source term in the transport equations, however ex- 
plicitely as the presence of well defined helicity bands. The transport equation takes 
the form 

iTiv%-V R g+leT3-A\g}o=6 (57) 

with normalization condition gog = —n 2 l. Here, the velocity is strongly renor- 
malized due to spin-orbit interaction. The quasiparticle trajectories are different for 
different helicity, and no coherence exists between the different helicity states. The 
matrix dimension can be reduced by a factor 2 compared to the case of weak spin- 
orbit splitting, and instead the number of Fermi surface sheets is increased by a 
factor of 2. If measurements are made that are spin-selective, the corresponding 
vector of Pauli spin matrices must be transformed according to 

Oil' = (UtOUfrw = U kXa a ap U^ ip . (58) 



3.4 Riccati parameterization 

One of the main obstacles of quasiclassical theory has been the non-linearity that 
is introduced by the normalization condition. A powerful way to deal with this 
problem is the choice of a parameter representation that ensures the normalization 
condition by definition. In this representation, the Keldysh quasiclassical Green's 
function is determined by six parameters in particle-hole space, y R A , y R,A ,x K , x K , 
of which •y R ' A , yR'A are the coherence functions, describing the coherence between 
particle-like and hole-like states, whereas x ,x are distribution functions, describ- 
ing the occupation of quasiparticle states [20 . 70 1 . The coherence functions are a 
generalization of the so-called Riccati amplitudes [70-74] to non-equilibrium situ- 
ations. All six parameters are matrix functions with the dimension determined by 
the degeneracy of the quasiparticle trajectories, and depend on Fermi momentum, 
position, energy, and time. The parameterization is simplified by the fact that, due to 
symmetry relations, only two functions of the six are independent. The particle-hole 
symmetry is expressed by the operation X which is defined for any function X of 
the phase space variables by 

Q(p F ,R,z,t) = Q(-p F ,R,-z*,t)*. (59) 

Here, z — £ is real for the Keldysh components and z is situated in the upper (lower) 
complex energy half plane for retarded (advanced) quantities. Furthermore, the sym- 
metry relations 

r A = (r R ) t , r A = (r R ) t , * K = (* K ) + (60) 
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hold. As a consequence, it suffices to determine fully the parameters y R and x K . 

The quasiclassical Green's function is related to these amplitudes in the following 
way [here the upper (lower) sign corresponds to retarded (advanced)]: 

5 r.a_^,._ ({i-y°?y l °{i + yoy) 2{\-yoy)- l oy \ R ' A 

8 + \ -2(l-foy)- l oy -(l-Yoy)- l o(l+yoy) J ' {01> 

which can be written in more compact form as ll82l 

/ eg & \ RA 
| R ' A = T 2ni f _ # A ±int 3 , (62) 

with the abbreviations r S = (1 — y°y)~ l and ^ = ^07. For the Keldysh component 
one can write |[82l 

Here, the o-symbol includes a time convolution as well as matrix multiplication; the 
inversion is defined with respect to the o-operation |82|. 

From the transport equation for the quasiclassical Green's functions one obtains 
a set of matrix equations of motion for the six parameters above [20 . 70 1 . For the co- 
herence amplitudes this leads to Riccati differential equations |f72l , hence the name 
Riccati parameterization. 



3.5 Transport equations 

The central equations that govern the transport phenomena have been derived in 
Ref. [20 70 1 . The transport equation for the coherence functions y(p F ,R,e,t) are 
given by 

(iSv F - V /f + 2e)y R ' A = [yoAoy+Eoy-yoE-A] R - A . (64) 

For the distribution functions x(p F ,R, e,t) the transport equations read 

(ihv F -V R + ihd,)x K - [yoA+Z] R ox K -x K o[Aoy-Z] A 

= -y R oE K oy A + A K oy A + y R oA K -E K . (65) 

The equations for the remaining components are obtained by the symmetry relation 
Eq.(E5J. 
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The transport equations must be complemented with boundary conditions for the 
coherence amplitudes and distribution functions at interfaces and surfaces 17511761 . 
For spin-active scattering such conditions were obtained in Ref. IT771 . Explicit for- 
mulations in terms of special parameterizations were given in Refs. 120, 781 4811 . 
Further developments include strongly spin-polarized systems [82-86], diffusive in- 
terface scattering lISTll or multi-band systems [88 1. We adopt the notation EOl that 
incoming amplitudes are denoted by small case letters and outgoing ones by capital 
case letters, see Fig. [8] Note that the velocity direction of trajectories is opposite 
for holelike and particlelike amplitudes as well as advanced and retarded ones. The 
boundary conditions express outgoing amplitudes as a function of incoming ones 
and as a function of the parameters of the normal-state scattering matrix. 




Fig. 8 Notation for the coherence amplitudes and distribution functions at an interface. Indices 1 
and 2 refer to the sides of the interface. The arrows for the Fermi momenta are for particle like exci- 
tations. The Fermi velocity directions are given by the directions perpendicular to the Fermi surface 
at the corresponding Fermi momentum. Quasiparticles move along the Fermi velocity directions 
(dashed lines). The components of the Fermi momenta parallel to the surface are conserved (indi- 
cated by the thin dotted line). For each trajectory, small case letters denote coherence functions and 
distribution functions with initial conditions from the bulk, and capital case letters denote functions 
with initial conditions at the interface. The interface boundary conditions must express all capital 
case quantities in terms of the small case quantities. Here, the simplest case, that involves only one 
Fermi surface sheet on either side ('two-trajectory scattering'), is shown. After Ref. 1201. 
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3.6.1 Coherence amplitudes 

The boundary conditions for the coherence amplitudes are formulated in terms of 
the solution of the equation ll82l 

fok] R = JX°#°^ (66) 
p 

n^f = [ie + 1 o % o fa] R (67) 

(the trajectory index p runs over all incoming trajectories) for [rj^-jf ] R , where the 
trajectory indices k,k! ,k\ run over outgoing trajectories involved in the interface 
scattering process, and the scattering matrix parameters enter only via the "elemen- 
tary scattering event" [}^] R,A - The quasiclassical coherence amplitude is given by 
the forward scattering contribution of [r k ^ k i] R , 

*i R =j£*- (68) 

Analogous equations [82 1 hold for the advanced and particle-hole conjugated com- 
ponents, [f p ^ pl ] R , [r pl ^ p ] A , and [f k/ ^ k ] A . 



3.6.2 Distribution functions 

For the Keldysh component not only the forward scattering contribution of [i^-f ] R 
is required, but also the off-scattering part 

[T k ^} R = {r k ^ k ,-r k 8 kk ,] R . (69) 

The boundary conditions for the distribution functions read 1 82 ] 

[^] K = E5 R ox K o5 A ,. (70) 

p 

x k K = E fe. + r ^-i ° %J R ° Kk 2 f ° [Shk + 7k 2 oT k2 ^k] A 

-E^j'oiMF^/, (71) 
h 

which depends on the scattering matrix parameters only via the elementary scatter- 
ing event [jcj^,] . Analogous relations hold forX^. 

The transport equation for the distribution function is solved by any function of 
energy in equilibrium. The correct boundary conditions in this case are 



= (1 
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for excitations of charge q in an electrostatic potential <P. 
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3.6.3 Case 1: one-trajectory scattering 

In this case only one incoming and one outgoing trajectory are coupled via the 
boundary conditions. The corresponding normal state scattering matrix is denoted 
by S and is a scalar in trajectory space. The boundary conditions read in this case 
simply 

[/] R = 5 R o 7 R o5 R , r R -[/] R , (73) 

and 

[x'] K = S R o/ o S A , Xf = [x'] K . (74) 

3.6.4 Case 2: two-trajectory scattering 

This is the case of scattering from two incoming trajectories into two outgoing tra- 
jectories. Examples are reflection and transmission at an interface, or reflection from 
a surface in a two-band system. The scattering matrix and the elementary scattering 
events have in this case the form 

V 21 22 / 1=1,2 1=1.2 

We give the solutions for trajectory 1, the remaining solution can be obtained by 
interchanging the indices 1 and 2. The boundary conditions read for i, j = 1,2 

^1=^1+^20^0)4^, ^=[^+^20^0/2/. (76) 

The equation for the ri«_2 can be solved by simple inversion, 

^2=W2°(l-fco^r 1 ] R , (77) 
and the solution introduced into the equation for r^_ l — Fj R , 

fi R = [Yn+Yi2°(l-?2oY 22 )- 1 or 2 oy 21 ] R . (78) 

For the distribution function one needs the components rf^ 2 = an d obtains 

Xf- = [*' n ] K + m 2 o f o [x 21 ] K + [x' 12 ] K o ^of*, 

+r£-2 ° (ft ° Mlf °tf-%)°f£,y. (79) 
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We present here formulas for the special case of the zero temperature conduc- 
tance when a single band system is contacted by a normal metal. We assign the index 
1 to the normal metal side of the interface and the index 2 to the superconducting 
side. The momentum for incoming trajectories on the superconducting side of the 
interface is denoted by k 2 , and that for the outgoing trajectory on the superconduct- 
ing side by k 2 . For the normal side the corresponding momenta are k\ and k\ (see 
Fig. [8] for the scattering geometry). The projection on the interface of all four mo- 
menta is equal. The corresponding incoming coherence functions in the supercon- 
ductor are 72(e) = )f(k2 ) e) and 72(e) = y 2 (k 2 ,e). Furthermore, Sn = Sf 2 (k u lc2), 
S22 = S 22 (k 2 ,k2), and 622 = ^22(^2,^2)- The Fermi velocity for outgoing direc- 
tions on the normal side will be denoted by Vfi = v>fi (k^). Having thus specified all 
momentum dependencies, we will suppress in the formulas below the momentum 
variables. In the case under consideration, after introducing Eqs. (l77l >. d78l . and d79l > 
into Eq. (1531 1. we obtain after some algebra (we omit hereafter the o sign) 



G(eV) 
Gn 



nv F i 

+ (nv F1 



Si2[l+A 2 (e)S 2 2] " S 12 A 2 (e) 
S 12 [l+A 2 (e)S22]72(£)52i 



2\ + 

£=-eV 



£=eV 



(80) 



where 



A 2 (e) = 1 - 72(£)S227>(£)S22 y2{e)S 2 2fi(e) 



-1 



(81) 



and we used the notation ||A|| 2 = |Tr(AA T ) for any 2x2 matrix A. The symbol 
(• • - )t=eV denotes Fermi surface average only over outgoing directions, and the ar- 
gument is to be taken at energy eV . For S22 = $22 = — y/R(d), S12 = S21 = \/D(9) 
(with impact angle 0), Eq. (f80T > reduces to Eq. ( fTob . For the tunneling limit we can 
neglect the second line in Eq. d80l >. and using the relation 



1+A 2 (e)5 22 = i{^(e) + 1} 



(82) 



with the complex quantity 



^i(e) = |(l -72(e)52272(e)52 2 ) ' (l + Ji^iU^Sn 



the conductance simplifies after some re-arrangements to 

= iReTr(nv F1 {S l2 ^2(eV)(S 12 ) f } 



(83) 



For the tunneling limit, in JV2 the surface scattering matrix (i.e. for S\2 = 52 1 = 0) 
can be used, for which the local density of states at an impenetrable surface is 



^ = iReTr(^ 2 (e) +5 2 2^2(e)4 

iV2 F 



(84) 
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3.6.5 Scattering matrix for non-centrosymmetric/normal-metal junction 

For the case that a non-centrosymmetric material with small spin-orbit splitting is 
brought in contact with a normal metal, we can use the formulas of the last sub- 
section. The scattering matrix for scattering between the two helicity bands in the 
non-centrosymmetric metal (index 2) and the two spin bands in the normal metal 
(index 1) can be expressed in terms of the scattering matrix for scattering between 
spin states on both sides of the interface. The corresponding transformation is 



S'n S'n\ = f 1 °\ ( s " s n\ f l 
S' 2l S' 22 J \pU k J'\S2iS22j'\pul 

s' 21 s' 22 J [pu*_ k )\s 2l s 22 J\ou^ 



(85) 
(86) 



For the simple case of a spin-conserving scattering in the spin/spin representation, 
the spin/helicity representation of the scattering matrix takes the form 





s[ 2 ) 




\s' 21 


$22 ) 


\t*U k 



tUl \ (~S'n~S\ 2 \_( r t*Ul k _ 



rUkUl)' \S' 2l S' 22 J \tUt k -rW_ k Ul k 



■(87) 



where r = and t = with r 2 + \t 2 = 1 are reflection and transmission coef- 
ficients that depend on the (conserved) momentum projection on the interface. We 
have chosen r real, as in quasiclassical approximation possible reflection phases do 
not affect the results. The case t — can be used to describe scattering at a surface. 

In the case of a contact with a non-centrosymmetric metal with strong spin-orbit 
split bands the scattering matrix has a more complicated structure. It connects in 
this case three incoming with three outgoing trajectories, and the scattering at the 
interface will not be spin-conserving. For this case, it does then not make sense 
anymore do use a spin/spin representation, but a spin/helicity representation must 
be used consistently. The scattering matrix must be obtained in agreement with the 
symmetry group of the interface, and it cannot in general be related anymore to the 
Ilk matrices in a simple way. 



3.7 Superconducting order parameter 

For the case of weak spin-orbit splitting one expects that to leading order in the 
small expansion parameters either a singlet or a triplet component nucleates. On 
the other hand, any finite spin-orbit interaction leads to a mixture of spin singlet 
(A s ) and triplet (A t ) components |89,90|. Consequently, the singlet or triplet states a 
never pure, but they are mixed. This mixing becomes in particular prominent when 
the spin-orbit interaction is strong. In this case, it does not make sense anymore to 
speak about singlet or triplet components, but it is necessary to start from the helicity 
basis. 
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It is interesting to consider what happens in the weak case first. In this case the 
triplet component is expected to be induced directly by the structure of the spin-orbit 
interaction, and the spin triplet component aligns with g(k). The gap function is in 
this case in spin representation given by, 

A° = (A k +D k g(k)-o)io® (88) 

which transforms in helicity basis into 

A = U k (A k +D k g{k)-a)ia^U T k 
= U k (A k +D k g(k) ■ a)U T k U k ia^U T k 
= (A k +D k \g(k)\a^)U k U^ 2 K (89) 

We introduce the notation 

(u k ul k ) u , = ( _° e 7 S ) - (90) 



Note that the identities (a^) 2 = 1, = -o®, and o^a^a® = -a^>, 

hold. With this, we can obtain the Nambu-Gor'kov space structure of the order pa- 
rameter 



A = 



0A\ = f {A k +D k \g\o 3 )a&aW 

A0)~ {(A_ k + D_ k \g\a^)*a^*a^ 

(A k + D k \g\a^)ia^\ f ia { ^ 

(A* k +D* k \g\G^)ioW )\ -ia(-*)* 



.(91) 



where A_ k = A k , and D k = D k , and we have used a^*a( 3 'cr^* = -a' 3 '. With 
A±(k) = A k ±D k \g\ the order parameter can be cast in the form 

with phase factors ^(Jfc) = -e,- a <Ps. Note that t x {-k) = -tx(k), and \t x {k)\ = 1, 
and A±(-k)=A±(k). 

We note that other possibilities to define the canonical transformation that diagonalizes the 
kinetic part of the Hamiltonian exist, which differ by the relation between particle and hole com- 
ponents. Using these alternative definitions (e.g. in Refs. |36 44 1), the order parameter is purely 
off-diagonal instead of diagonal in the band representation, and the symmetry relation Eq. ([59) 
becomes non-trivial (see e.g. Ref. |44|). Here, we prefer a transformation that preserves the sym- 
metry \59\ , and renders the order parameter above diagonal. This is a natural choice when treating 
strongly spin-orbit split systems, where the order parameter should be band diagonal. 
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The coherence amplitudes in a bulk system with order parameter Eq. d92l are of 
a similar form, 

*M^+n ),+< \_ (t ,°),_(*)) 
KM, - («n' +Mr w M )°M- t) .) 

with f±(k,e) = Y±(—k,—e)*. In inhomogeneous systems helicity-mixing can take 
place. If this happens, the form of the coherence functions is the same band-diagonal 
form as above for the case of strong spin-orbit splitting, however has the full matrix 
structure for the case of weak spin-orbit splitting. 



3.8 Results 



3.8.1 Andreev bound states near the surface 

The surface bound states are determined by the poles of the Green's function. Fol- 
lowing Refs. [42 , 44], we consider specular reflection, whereby the component of 
k normal to surface changes sign, k — > k, whereas the component parallel to the 
surface is conserved. We find the amplitudes y(k, e) by integrating forward along 
the incoming, k, trajectory starting from the values in the bulk, and the amplitudes 
f(k, e) by integrating backward along the outgoing, k, trajectory, again starting from 
the values in the bulk [20 1 . For the homogeneous solutions one obtains 



iyJ\A±{k)\ 2 -e 2 £ + i^\A ± (k)\ 2 -e 2 



Note that the spin-orbit interaction in the helicity basis enters as a term proportional to c' 3 ', 
see Eq. <55K Consequently, this term commutes with any term diagonal in the helicity basis, and 
thus drops out of the homogeneous solutions in Eq. \96\ (see Ref. ( 9Tj for the case of a Rashba- 
type spin-orbit coupling). Note, however, that this is not in general the case for non-homogeneous 
solutions: when helicity mixing takes place due to impurities or surfaces and interfaces, and a fully 
self-consistent solution is obtained, then the spin-orbit coupling term in Eq. {55j enters through the 
transport equation. 

The amplitudes fjt and f^, are determined from the boundary conditions at the 
surface. We consider here a simple model of a non-magnetic surface, that conserves 
the spin under reflection (this assumption only holds for a small spin-orbit interac- 
tion in the bulk material). In this case the components of g in the spin basis, 

g°(k,e) = OZg(k,e)t) k , (97) 



are continuous at the surface. This leads to a surface induced mixing of the helicity 
bands according to 
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ulr%e)U*_ k = r(k,e) = f(k,e) = ulr(k,e)Ul k , (98) 
(k,e)U k =f s (fee) = f (fee) = l£ s y(fce)CV (99) 

Note that these boundary conditions correspond to Eq. (F73l > with 5 R and 5 R given 
by the (22)-components of Eq. d87] >. 

We proceed with discussing the local density of states at the surface, N(e), that 
is defined in terms of the momentum resolved density of states, N(k, e) by 

N(k,e)/N F ^-(27z)- 1 lmTr l {g(k,e)}, N(e) = (iV(fee)}, (100) 

which can be expressed in terms of the coherence amplitudes in the following way 
(here k points towards the surface and k away from it), 

N(k,e)/N F = ReTa { [l - y(fe e)f (fee)] _1 - 1/2} 

^V(fee)/^V F =ReTr^{[l-r(fee)7(fee)]- 1 -l/2}. (101) 

We obtain r(k, e) and f (fee) from Eqs. d98l l and (|99l , with /(fee) and f(fe e) from 
Eqs. (HO), 495J, @, and ([T05l 

The bound states in the surface density of states correspond to the zero eigenvalues of the matrix 

1 - r(k,e)f (fee) = 1 - r(k,e)(Ut k U T t mk,e)(UkUk) ( 102 ) 

at the surface. An explicite calculation results in an equation for the Andreev bound states energy 
in terms of the surface coherence amplitudes in the helicity basis |44|, 

{l + y + y + }{l + y_y_} 

{i + r + r-}{i + y-r + } ' ( } 

where we used the abbreviations y± = y±(k,e) and y± = y±(t, e). The "mixing" factor ^ is 
determined by the change of g(k) — > under reflection k — > k at the surface, 



n 2^_a + sin 2^a t an 2 



- cos z -V- 2, tan- 



(104) 



where 9 g , (p g and g ,(p g are the polar and azimuthal angles of g(k) and g(/fc), respectively. 



In general, the order parameter must be obtained self-consistently at the surface. 
Helicity mixing at the surface will lead necessarily to a suppression of the order 
parameter. To gain insight in the role of the order parameter suppression it is useful 
to model it by a normal layer of width W next to the interface. Trajectories incident 
at an angle a* from the surface normal travel through a normal region of an effective 
width 2W* = 2W/cos(o!jt). Thus, the surface coherence amplitudes gain a phase 
factor, 

y±(fee) = f ± (k,e)e 2i£W ^ cos ^\ y±(fee) = f ± (k,e)e 2i£W ^ C0S ^ . (105) 
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Fig. 9 (Color online) Local surface density of states N(e) /Np for a Rashba superconductor, g(k) = 
cCrIc x z. The surface is parallel toz. The curves are for A± = Ao(q± \g(k)\)/(q + 1), with q ranging 
from to 2. In (a) and (c) the order parameter is assumed constant up to the surface, and in (b) 
and (d) a suppression of the order parameter to zero in a surface layer of thickness W = 2^o with 
<i;o = livf /27tks T c is assumed, (a) and (b) is for a cylindrical Fermi surface, vp = (v v , v y , 0) , and (c) 
and (d) is for a spherical Fermi surface. (The symbols are labels for the curves only). 



Similarly like for g we will use in the following polar and azimuthal angles for 
the vector k, defined by {k x ,k y ,k z } = |^|{sin(0 t )cos(%), sm (0t) sm ( ( Pt), cos (0Jt)} 
(where < < n). We also introduce the notation g = g/ma\(\g\). For the order 
parameter, we assume isotropic A^ = A and = D in Eq. d88l ). and introduce the 
parameter q = A /D' where D' = D ■ max(|g|) [42]. In this case A± — D'(q ± |£|) 
with maximal gap amplitudes Aq — D'(q+ 1). 

In Fig.|9]we show results for a Rashba-type spin-orbit coupling, g(k) = auk x z. 
The surface is aligned with the £ direction. In (a) and (b) we use a cylindrical Fermi 
surface, for which \g(k) \ = 1 . In (c) and (d) the results for a spherical Fermi surface 
are shown, for which \g(k) \ = sin(fy). The effect of a surface layer with suppressed 
order parameter is illustrated in Fig.|9](b) and (d), where Eq. (|105t with W = 2^q is 
used, where <^o is the coherence length = hv^/2nk^T c . 

For the special case q = we have A + = — 4_ = 4osin(0^) (we use a real gauge). For this 
case. 0g = Of, = n/2, and consequently, . // = tan 2 (p g = cot 2 <£><-. The bound states are then given 

bv B2ll44l 

e ( 2W £ \ 

— = -sin -±<p k sin(e t ). (106) 
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Numerical solution of the problem shows that the "principal" bound state branches 
e((pic) with energies away from the continuum edge contribute the most to the sub- 
gap DOS. For W ^ the main branch develops a maximum at £* < Aq, 
which gives rise to a peak in the surface DOS near £*, see Fig.|9](b) and (d). Fully 
self-consistent solution confirms this fl4l . For q — > °° the order parameter becomes 
insensitive to helicity mixing, i.e. the effective W decreases for increasing q. 

Andreev bound states in non-centrosymmetric superconductors have unusual 
spin structure [40 44 |. It is found, that the states corresponding to different branches 
of Eq. (11061 > have opposite spin polarization. Since the spin polarization changes 
sign for reversed trajectories, the Andreev states carry spin current along the inter- 
face. Such spin currents exist in NCS materials because the spin is not conserved, 
and consequently precession terms enter the continuity equation 1 92 1 . There are spin 
currents both in the normal state and in the superconducting state. As was found in 
Ref. PBl . the most prominent feature is a large surface current with out of plane spin 
polarization (reminiscent to that in spin Hall bars [93 |) that flows along the surface, 
and decays rapidly into the bulk on a Fermi wavelength scale. In addition, there is 
also a surface induced superconducting spin current with out of plane spin polar- 
ization, that adds to the background microscopic spin currents and greatly exceeds 
them in the limit of small spin-orbit band splitting. This effect is in this case solely 
determined by the structure of the superconducting gap. Superconducting spin cur- 
rents decay into the bulk on the scale of the coherence length and show oscillations 
determined by the spin-orbit strength due to Faraday-like rotations of the spin co- 
herence functions along quasiparticle trajectories [44|. 

3.8.2 Tunneling conductance 

For a three-dimensional model, which for the Rashba-type spin-orbit coupling was 
discussed in Ref. [42], we present in the following tunneling conductances in various 
geometries. We will discuss several types of spin-orbit interaction: 



For the symmetry C^ v , corresponding to the tetragonal point group, the two pa- 
rameters 77 and tj' can both be non-zero. We will discuss below the special cases 
77 = and 77/ = 0. The case tj' = corresponds to a Rashba spin-orbit coupling. 
The type of spin-orbit coupling we consider for the full tetrahedral point group, Tj, 
is also known as Dresselhaus coupling. Finally, for the cubic point group, 0, the 
simples form for g is considered here, which is fully isotropic. All the cases above 




6 


k x kyk^{k x 



: 
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are relevant for non-centrosymmetric superconductors: C$ v for CePt3Si, CeRhSi3, 
and CelrSi 3 , Id f° r Y2C3 and possibly KOs 2 6 , and for Li2(Pdi_^Pt v )3B. 

The zero temperature tunneling conductance is obtained according to the formula 
Eq. d83l l. which leads for a spin-inactive 5-function barrier to 



where a k is the angle between the surface normal and k. A remark is in place here. In 
principle, the interface barrier will be spin-dependent once a spin-orbit split material 
is brought in contact with a normal metal. However, for the limit of small spin-orbit 
splitting we can neglect the spin-dependence of the interface potential consistent 
with the quasiclassical approximation. The corrections are of the same order as the 
corrections for the quasiparticle velocity in this case, and are of higher order in the 
parameter SMALL. 

In Figure. [10] we show the tunneling conductance G(eV)/G^ obtained from 
Eq. ( |108b with Eqs. J lOOt - dToTb for various types of spin-orbit coupling corre- 
sponding to the spin-orbit couplings in Eq. ( 1107b . and for various alignments of 
the surface normal with respect to the crystal symmetry directions. The tunnel- 
ing parameter in this figure is Do = 0.1. We show curves for an order parameter 
A± =Ao(q± \g(k) )/{q + 1 ), with q ranging from to 2. For simplicity, we concen- 
trate here on the assumption that the order parameter is constant up to the surface, 
i.e. we use the bulk solutions Eqs. (|94b-(l96l). For a detailed quantitative description, 
a self-consistent determination of the order parameter suppression near the surface 
must be obtained. We also use the simplifying assumptions of a spherical Fermi 
surface with isotropic Fermi velocity. 

As seen from Fig. [TO] a rich structure of Andreev bound states below the bulk 
gap energy develops, that depends strongly on the alignment of the surface with 
the crystal symmetry axes. In (a) and (b) a pure Rashba spin-orbit coupling lc = 
[ky,—k X) 0] on a Fermi sphere is assumed. Below the critical value q = 1, a zero 
bias peak appears for tunneling in the direction perpendicular to the z direction, i.e. 
the (1,0,0) or (0,1,0) direction, however a dependence quadratic in energy appears 
for tunneling parallel to the z direction, i.e. the (0,0,1) direction |42|. For q > 1 the 
tunneling density of states acquires a gap, as then the singlet character of the order 
parameter dominates. 

In Fig. [TO] (b) and (c), we show results for a hypothetical spin-orbit coupling of 
the form g = rj'[0,0,k x k y k z (k^ —fey)], that is consistent with the same point group 
symmetry D4 V as the Rashba spin-orbit coupling. For this case, a sharp zero bias 
conductance peak exists for q < 1 in all tunneling directions. In contrast, for q > 1, 
the zero bias conductance peak only exists when tunneling perpendicular to the z- 
direction, however, not when tunneling parallel to the z-direction. 

In Fig.[l0](d) we consider the cubic point group symmetry, and assume the sim- 
plest form of a fully isotropic spin-orbit interaction of the form g = \k x ,k y ,k z ] . Here, 
for q < 1 the tunneling conductance is zero at zero bias, but raises sharply away from 



G(eV) _ (cos(a k )D(a k )N(k,eV)) 
G N ~ (cos(a k ) D{a k )) 



D{a k ) 



Do cos 2 (a*) 
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(a) C 4v (n'=0) , 
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Fig. 10 (Color online) Tunneling conductance C(eV) /Gn for various types of spin-orbit coupling 
corresponding to the indicated symmetry groups, and for various alignments of the surface normal 
h as indicated. The spin-orbit vector is of the form C^ v : g = T] [k y , —k x ,0] + r]'[0,Q,k x k y k z (k 2 — k^)]\ 

0:g = [h,k y ,k z ] (this case is fully isotropic); T d : g = 2[k x (k 2 - B),k y (& -k 2 x )X(il -£?)]• The 
curves are for A± = Ao(q ± \g(k)\) / {q + 1), with q ranging from to 2. The order parameter is 
assumed constant up to the surface, and a spherical Fermi surface with isotropic Fermi velocity is 
assumed. The tunneling parameter is Do = 0.1. Curves are vertically shifted by multiples of 0.2. 



zero bias, showing side peaks due to Andreev bound states. At q = 1 this structure 
disappears with only a pseudogap remaining. For q > 1 a gap opens. 
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Finally, in Fig. [10] (e)-(g) we show results for the full tetrahedral point group 
Trf. We compare tunneling in (1,0,0), (1,1,0), and (1,1,1) directions. Note that in 
this case, the relation (11071 ) between g and k is not invariant under a rotation of 
both vectors by 90 degree around the k x -, k y -, or L-axis, but an overall sign change 
appears; however, the conductance spectra are insensitive to this sign change. For 
q < 1 there is a vanishing zero bias conductance for tunneling in (1,0,0) direction, 
and a low-energy dispersive Andreev bound state branch for tunneling in (1,1,0) 
direction. For tunneling in (1,1,1) direction, the zero bias conductance vanishes for 
q = 0, and shows a sharp zero bias peak for < q < 1 . For q > 1 the tunneling 
conductance becomes gapped for all directions. 

As can be seen from these results, studying directional resolved tunneling in non- 
centrosymmetric superconductors gives important clues about the order parameter 
symmetry and the type of spin-orbit interaction. 

3.8.3 Andreev point contact spectra 

Here we present results for the case of a point contact between a normal metal and 
a non-centrosymmetric superconductor. We use Eq. ( f80l > to calculate the spectra, 
with a scattering matrix that has the form shown in Eq. (187) . We assume isotropic 
Fermi surfaces in the materials on both sides of the interface, and for simplicity use 
equal magnitudes for Fermi momenta and velocities. The transmission amplitude is 
modeled by that for a 5-function barrier, 



and the component of the Fermi velocity along the interface normal in direction of 
current transport is hvpi = v'FCOs(at). 

In Fig.Qj] the Andreev conductance G(eV)/G-^ for various types of spin-orbit 
coupling and for various alignments of the surface normal h are shown. Here, the 
transmission probability Dq = ?q is varied from zero to one. We restrict here to the 
case q = 0, i.e. an order parameter of the form A± = ±Ao\g(k)\). Again, a spherical 
Fermi surface with isotropic Fermi velocity is assumed. We also compare the case of 
a surface layer with suppressed order parameter (dashed lines) with that of an order 
parameter constant constant up to the surface (full lines). To model the order param- 
eter suppression, we assume a layer of thickness W = 2E,q with E,q = hvp j2%k^T c 
(dotted lines) in which the order parameter vanishes. Thus, we use Eq. ( 1105b as 
incoming solutions for the coherence amplitudes. Note that for Do = 1 the surface 
layer with zero order parameter does not affect the Andreev conductance. This is due 
to the fact that for perfect transmission the normal region simply extends slightly 
further towards the superconductor, and within our approximation we neglect the 
spin-orbit effects in interface potential. 

For a larger spin-orbit coupling the interface between a normal metal and a normal conducting 
non-centrosymmetric metal with strong spin-orbit interaction becomes necessarily spin-active, as 



f COs(0Ejt) 



(109) 
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Fig. 11 (Color online) Andreev conductance G(eV)/Gff for various types of spin-orbit coupling 
corresponding to the indicated symmetry groups, and for the indicated alignments of the surface 
normal h. The transmission probability Do = is varied. The spin-orbit vector is of the form 
Qv : S — nfei —k x ,0] + i'[0,0,^jLkj(^ — £?)]; O: g — [kx,ky,k z ] (this case is fully isotropic); 
T d : g = 2[k x (i 2 -k 2 ),k y (k 2 - k 2 x ),k z (il -%)]• The curves are for A± = ±4 |^(fc)|). The order 
parameter is assumed constant up to the surface (full lines) or suppressed to zero in a surface 
layer of thickness W = 2ijo with £0 = ftvp /27tkg T c (dotted lines). For Do = 1 these two cases give 
identical results. A spherical Fermi surface with isotropic Fermi velocity is assumed. 



the interface potential term in the Hamiltonian must be hermitian. Thus, a perfect transmission is 
not realistic in such a case. For weak spin-orbit splitting these effects are also present, however 
modify the results only to order vso/£f> or on energy scales v| /£p. 
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For lower transmission, we remark as an overall observation that the suppression 
of the order parameter does not affect the value of the Andreev conductance at zero 
bias. This is simply due to the fact that in the clean limit the coherence amplitudes 
become effectively spatially constant for e = 0. For higher bias, deviations can be 
observed, that in general lead to a shift of Andreev bound states to lower bias. 

We turn now to the Andreev point contact spectra for frj = 1 ■ As can be seen, the 
form of the spectrum is sensitive to the type of spin orbit coupling, and the associ- 
ated order parameter symmetry. For a Rashba spin-orbit coupling, Fig.[TT](a) and 
(b), the Andreev conductance is enhanced to twice the normal conductance at zero 
bias, however to a smaller value for finite bias. There is a pronounced anisotropy 
in the shape of the Andreev conductance spectra. In (c) and (d) the Andreev con- 
ductance shows a sharp kink feature at zero bias, associated with the complex nodal 
structure of the spin-orbit vector. For cubic symmetry, (e), we observe an Andreev 
conductance resembling that of an s-wave spin singlet superconductor. And, finally, 
for a Dresselhaus spin-orbit coupling, (f)-(h), the Andreev conductance shows a be- 
havior similar to the case of a Rashba spin-orbit interaction, however with a not so 
pronounced anisotropy. 



4 Conclusions 

We have given an overview over the current status of the theoretical understanding 
of Andreev bound states at the surface of a non-centrosymmetric superconducting 
material, and have presented results for tunneling conductance, point contact spec- 
tra, and spin polarized Andreev bound state spectra. 

The new feature in non-centrosymmetric superconductors is the possible appear- 
ance of spin polarized Andreev states, that carry a spin-current along the interface or 
surface. The presence of such Andreev bound states that cross the chemical potential 
as a function of incident angle to the surface, is a topologically stable superconduct- 
ing property. Such bound states exist as long as triplet order parameter components 
(in spin representation) dominate singlet components of the order parameter. When 
both components are equal, the bound states at the chemical potential disappears, 
and a topologially new ground state establishes. The transition between the two 
states is a quantum phase transision. 

The spectrum of Andreev states at the surface provides valuable information 
about both the structure of the superconducting order parameter and the vector field 
of spin-orbit vectors on the Fermi surface. In this chapter we have concentrated on 
the rich structure that appears for the limiting case of a small spin-orbit splitting of 
the energy bands in the non-centrosymmetric material. In this limit, the spin quan- 
tum number is approximately conserved during scattering from surfaces and inter- 
faces with normal metals, which leads to strong mixing between the helicity bands 
in the non-centrosymmetric material. The opposite limit of strong spin-orbit split- 
ting is still largely unexplored. We have provided a theoretical basis in this chapter 
that allows to treat this case as well. 
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Finally, we would like to mention that interesting effects, like e.g. effects related 
to the spin Hall effect, or to Berry phases associated with the change of the spin-orbit 
vector along closed paths, are interesting subjects left for future studies. 
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